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THE OBERBECK BOUSSINESQ APPROXIMATION 
AS A CONSTITUTIVE LIMIT 

YOSHIYUKI KAGEI AND MICHAEL RUZICKA 


Abstract: We derive the usual Oberbeck-Boussinesq approximation as a con¬ 
stitutive limit of the full system describing the motion of an compressible lin¬ 
early viscous fluid. To this end the starting system is written, using the Gibbs 
free energy, in the variables v, 6 and p. The Oberbeck-Boussinesq system 
is then obtained as the thermal expansion coefficient a and the isothermal 
compressibility coefficient (3 tend to zero. 


1. Introduction 

The well-known Oberbeck-Boussinesq m, la approximation was designed as a 
simplified model for the thermo-mechanical response of linear viscous fluids under¬ 
going isochoric motions in isothermal processes but not necessarily isochoric ones 
in non-isothermal processes. Its roots stem from the end of the 19*^ century. Nev¬ 
ertheless, its justification from the point of view of continuum mechanics was quite 
recently given in 1996 by m (cf. m, 0 and m for earlier contributions). 

From the mathematical point of view the expansion used in m is still formal. 
In [5], [5] a rigorous justification of simplified problems has been given. We refer 
the reader to for a completely different approach in which singular limits of 

the full system are discussed. 

In this paper we use a new approach. Motivated by the studies in n, 0 we 
obtain the Oberbeck-Boussinesq approximation as a constitutive limit. Since this 
limit is singular it depends on the way how it is achieved. To achieve our result 
we re-write the full thermo-mechanical system with the help of the Gibbs free 
energy in the variables v, 0 and p. It is important that the Gibbs free energy 
depends on two (dimensional) parameters a, b, which tend to zero. Then we non- 
dimensionalize the system and make assumptions on the behaviour of the thermal 
expansion coefficient, the isothermal compressibility coefficient, the specific heat 
coefficient at constant pressure and the ’’density” on the non-dimensional versions 
of the paramaters a, b, which are denoted by A, B. We show that these requirements 
can be fulfilled by an easy example, which assumes that the density is linear in the 
pressure and the temperatur. Provided that weak solutions of the full thermo¬ 
mechanical system satisfying a uniform estimate exists we show that the limits 
satisfy the Oberbeck-Boussinesq approximation. Moreover, if we assume that the 
approximation parameters are fixed and sufficiently small we can also recover the 
results obtained by a power series expansion in |I2j . 

Let us finally introduce some notation: 
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In what follows boldfaced minuscules always stand for vectors and vector valued 
functions whereas boldfaced capital letters represent tensor valued functions, i.e. 
V = {vi,V2,V3y , T = := All quantities are considered 

at points x = (xi,X 2 ,xa)"'" S and at a certain time t. We use the abbreviations 
dk ■= dt ■■= (analogously for dg...),d'^ := dgdg, div v ;= diVi+d2V2+dzV'i 
and Vp := (i9ip, 92P, )^. The dot between two quantities denotes the corre¬ 

sponding scalar product, whereas the superposed dot is the usual material time 
derivative: 


V • w 


T L 


3 

Y, VkWk, 0 

k=l 

Y TkiLki, b 

k,l=l 


5(0 -I- V • V6» = dtO + Y VkdkO, 

k=l 

dth + [Vb]v = citb -f Y 'i’kdkh. 

k=l 


The trace of some tensor D is denoted by tr D and ||D|p :=D-D. For the identity 
tensor we write I. 


2. Derivation of the Approximation 

2.1. Governing Equations and Assumptions. The starting point for our anal¬ 
ysis is the balance of mass, linear momentum and energy and the Second Law of 
Thermodynamics in the form of the Clausius-Duhem inequality: 


p + p div V 

= 0, 


pv - 

divT 

= pb. 

(2.1) 

pe - 

- divq 

= T • D -1- p r. 


(!) 

r 

-P 9 

IV 

0 

(2.2) 


where p denotes the density, v the velocity field, T the symmetric Cauchy stress 
tensor, b the density of external body forces, e the specific internal energy, L the 
velocity gradient, r the radiant heating, 9 the temperature, rj the entropy and q 
the heat flux vector. 

In the following we neglect radiant heating, i.e. r = 0, and assume that the body 
forces has a potential, i.e. b = V/. Moreover, we restrict ourselves to the case of a 
compressible linearly viscous fluid. Thus, we assume that 


T = —pi -I- A (trD)I -t- 2pD, 
q = K V0, 


(2.3) 


where p is the pressure. A, p are the constant viscosities and k is the constant 
thermal conductivity. 

If V, 9 and p are considered as independent variables in (j2.1ll - ()2.3ll it is useful 
to introduce the Helmholtz free energy tf) through 

i;{p,9) ■=e{p,T])-9T]. (2.4) 

In this case we obtain from (12.31) that 

T] = -dgip , P = dpij) , 

p>0, 3A-l-2p>0, K>0, 


(2.5) 
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while ( 12 . 11 ) reads 


p + p div V = 0 , 

p \ — 2p, div D — A V(tr D) + V = P V/ , ^2 6) 

-p6({dli)) 6 + {dlsi)) p) - K A0 = 2/r |D| V A |tr D|" . 

However, for our purposes it is more convenient to view v, 6 and p as independent 
variables in To this end we introduce the Gibbs free energy (f) through 

(j){p,0) :=il;{p,e)+pp~K (2.7) 

In this situation we conclude from (j2.3|l that 


77 = -de(j), p ^ = dpcj), 
p>0, 3A + 2p>0, K>0. 


( 2 . 8 ) 


The system m now reads 

{dp(j))~^ (^{dgp(j)) 9 + (5p^)p) = divv , 

(ap()>)"V-2pdivD-AV(trD) + Vp= {dpct>)~Wf, (2.9) 

-9 {dp^y^ ((a|(^) 9 + {dly)p) -kA9 = 2p |D|" + a |trD^ 

Equation ( 1^ 1 nicely reflects the fact that changes in volume are induced by 
changes in temperature and changes in pressure. In fact, if we introduce the thermal 
expansion coefhcient a and the isothermal compressibility coefficient P through 

a{p,9) ■= {dp4>y^dlp(t), 

I3{p,0) := -{dp(l)y^d^(l), 

we can re-write (EH) 1 as 

aO — Pp = div V . 


( 2 . 10 ) 

( 2 . 11 ) 


It is well known that for many fluids the thermal expansion coefficient a is small 
(a Ri (10“^ — 10“^)K“^) and the isothermal compressibility coefficient /3 is even 
smaller (/3 « (10“^^ —10“^°)Pa~^). Finally, it is convenient to introduce the specific 
heat coefficient at constant pressure through 

Cpy9)-.= -9dl^. ( 2 . 12 ) 

With this notation we can re-write (EHa as 

Cp {dp(f)~^ 0 - a 6 lp - K A 6 » = 2/r |D| V A |tr D|^ (2.13) 

Remark 2.14. Note that in [12] special fluids that can sustain isochoric motions in 
isothermal processes have been considered. This formally corresponds to neglecting 
/3 in (12.111) . For such fluids the Oberbeck-Boussinesq approximation has been 
formally derived in |12| with the help of a power series expansion. In 0 , El a 
rigorous mathematical justification of a simplihed modell has been carried out. 
However, the mathematical justihcation starting with the full system from [12j is 
still lacking. One of the difficulties is the lack of appropriate apriori estimates, 
which is related to the fact that /3 has been neglected in (12.111) . 
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III [I] , [2] thermal expansion models have been considered as a constitutive limit 
for free energies. In these papers the predictions of the compressible theory has been 
compared to the prediction of different limiting theories. However, these limiting 
theories are from the mathematical point of view singular limits. Thus, the way 
how the limit is achieved is important and different ways can result in different 
limiting systems. 

In the present paper we combine ideas from m, 0 and m. m in order to 
derive the Oberbeck-Boussinesq approximation as a constitutive limit. To this 
end we consider the system (12.9|) for a family of Gibbs free energies depending on 
(dimensional) parameters a, 6 > 0, i.e. 

4'{p,9) = . 


The dependence on the parameters a, b will be suppressed in the notation in most 
cases. Under certain assumptions we consider the limit as a, b tend to zero. Be¬ 
fore that we render the system (12.91) non-dimensional by introducing dimensionless 
variables 


X := 


X 



V 



P 

5 

TT 



r 5 


/ 


gL^ 


where L,T,V,'!t and 0r,'d are typical length, time, velocity, pressure and tempera¬ 
tures, while g is the gravitational constant^ Moreover, let (po be a typical value for 
the Gibbs free energy. We will not choose V, L and T independent of one another 
but use the relation V T = L. Note that the constants /r, A, k, 6 r and g are assumed 
to be independent of the parameters a, b, while V, T, L, d, tt and po may depend on 
the parameters a, b. We will make this dependence explicit if necessary. 

The definition of the non-dimensional analogues of a and b is motivated by the 
following observation: We introduce non-dimensional a thermal expansion coeffi¬ 
cient a and the isothermal compressibility coefficient (3 by 


a{p, 0) 


a{p, d) 

ao 


/3{p, 9) 


13{p, 9) 
f3o 


where ao and /?o are typical values for for the thermal expansion coefficient a and 
the isothermal compressibility coefficient The non-dimensional version of (12.111) 
then reads 


ao'd a9 — PoT^ P P = div v, 


(2.15) 


where the superposed dot now stands for the non-dimensional material time deriv¬ 
ative. In typical applications we see that 

aoi? ~ 10”^ - 10"^ , Pott ss 10 "® - 10“® . (2.16) 


Since api? and /Sqtt are non-dimensional numbers (depending on the parameters 
a, b) we roughly want that the non-dimensional analogue A of a behaves as ao 'd 
and the non-dimensional analogue B of b behaves as Po tt. This is made precise in 


^ This form of the non-dimensionalization is motivated by a typical situation when the 
Oberbeck—Boussinesq approximation is used, namely a fluid layer in a gravitational field with 
a certain difference between the temperatur 9* at the top and 9^ at the botton. In this case we 
would choose 9 = B*" - 9^ and Or = -|- 9^). 
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the following way: We set 


h’^'{p,0) := 




±a,h 

K 


(2.17) 


and define the non-dimensional parameters A, B through 

A = Aia, b) := 1 ), 

B = B{a, b) := -(^^“■"(1, ^“•'’(1,1). 

Moreover, we assume that the mapping (a, b) i—^ (A, B) is invertible and that 
(a,b) —> (0,0) implies (A,B) —>■ (0,0). Finally, set 


(2.18) 


^ - MA,B),b(AB) 


(2.19) 


The system (12.91) for the non-dimensional quantities and differential operators then 
becomes (we skip all bars for convenience): 




{dpfj)^'^) ^V - (2p div D + A V(tr D)) + -^ Vp 

= ( 2 - 20 ) 


-{e + Or) {dp4>^’^) ^ 0 + A6» 

= ^(2/r|D|VA|trD|^), 

where we used the notation 0o := Of course, for the quantities 

V,T,L,'d and tt the same convention is used, e.g. V = yAA,B),b{A,B)^ 

Apriori there is no obvious representative velocity V in natural convection pro¬ 
cesses. As already observed in [4] such processes are reflected by the assumption 
fa g L ao This is translated in our situation by 

V^:=AgL. (2.21) 


From p. 201) 9 follows that non-trivial body forces are only possible if (j)o fa g L. 
Consequently, we set 


7 := 



( 2 . 22 ) 


and require that 7 = 0(1) as A, B tend to zero. We also introduce the Reynolds 
numbers Re^ and ReA as well as the Prandtl number Pr by setting 


Re^, := 


TT V L 


ReA := 


71 V J-J 


Pr := 




(2.23) 


2 p (/)o ’ ' A (/)o ’ 

and require that Re^^ = 0(1), ReA = 0(1), Pr = 0(1) as A,B tend to zero. Note 
that all these requirements can be fulfilled simultaneously, e.g. if 


V fa Ay Lfa A 3 


(j)o fa A 3 


TT 


A- 


■dfaA^ 


(2.24) 


then all above requirements are satisfied. 
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Using the above notation we can re-write (I2.20|) as 

= divv, 

(i?,0-v-^divD-^V(trD)-,5vp 
= (2.25) 


= A 


{e + Or) e + 

1 


1 


Pr Re.j 


■ Ae 


|D| 


|trD| 


^7Re;_, 7 Re A 

Concerning the behaviour with respect to A, B of the remaining quantities ’’thermal 
expansion coefficient”, ’’isothermal compressibility coefficient”, ’’specihc heat at 
constant pressure” and ’’density” in (12.251) we make the following assumptions: 
There exists constants co > i = 1,2, and functions /3j ’^{p,6), 

cf"^{p,6), pf’^{p,9), i = 1,2, such that 

dep(t>-^’^{p,0) 


dp(j)AB(^p^9) 

d^(j)^’^{p,9) 

dp(j)AB(j)^ 0 ) 

1 

dp<l)AB{j)^0) 


= A{1 + Aa^'^{p,9) + Ba2’^{j),9)) , 
= -B{l + Apt’''{p.0)+Bp^'^{p,9)) 


= (l_A{9 + 9r) + Bp + p^'^{p,9))k( 


A,B 


A,B 


-{9 + 9r)d^gct>^-^{p,9) = {co + A^cf’^{p,9) + B^ 4’^{p,9))k^' 
where we require that 

lim kf'^ = 1, 

(A.B)^(O.O) 

and locally uniformly in p, 9 

af'^{p,9)=0{l), /3f’''(p,0) = O(l), cf'^{p,9)=0{l), 

pf'^{p,9) = 0{A^)+0{B^) + 0{AB) 


(2.26) 


(2.27) 


(2.28) 


with * = 1, 2, as A, B tend to zero. 

Let us illustrate the above procedure and assumptions by the following example: 


Example 2.29. We assume that 

p°-'^(j),9) = {dp<j3°”’’{p,9)) ^ = po{l +bp-a9) , 

i.e. the density depends linearly on the pressure and the temperature. This consti¬ 
tutive equation is in dimensional form and the constants po, a, b are assumed to be 
positive. A possible Gibbs free energy compatible with such a behaviour is given 

by 

())“’^(p, 9) = b~^pQ^ (ln(l + bp — a9) — ln(l — a 0)) — cq9 (In 0 — 1), 
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where co > 0. From the definition in (12.1811 we easily compute 
A = A{a,b) = a'd{l-\-h-K — a'd{l + Or)) 

B = B{a^h) = hiril-\-hTT — ad{\ + Or)) 

We set 

1 ? := do (a do) ~ , tt := ttq (a do) ~ , (pg’’’ := </>“’*' (tt, , 

where do and ttq are positive constants independent of a, b. Moreover, we require 
that 

(2.32) 


(2.30) 


(2.31) 


lim -- -^77 

(a.&)-)-(o.o) (ado)^/^ 

Note that in typical situations this requirement is fulfilled (cf. (12.16^ 1. Moreover, 
this requirement ensures that {A,B) —>■ (0,0) if (a, 6 ) —>■ (0,0). The system of 
equations (12.301) for a, b is solvable and we compute 


= 0 . 


1 

a = — 


A 


\ 4/3 

do V ^ ^(1 T ^r) ~ B ) 


do \XA,B 


4/3 


1 

TTq 3 , 4/3 
^ ^A.B 


(2.33) 


Moreover, straightforward manipulations show (with skipped bars) 
dlp4>^'^{p,0) _ B{p-l)-A{0-1) 


dpip^’^ip 

d^cj)^’^{p,0) 

dp(j)^^^{p,0) 

1 

dp(t)^'^{p,0) 


^ ^ / B{p-1)-A{0-1) X 

i6») V i + B{p-l)-A{0-l))’ 


= -b(i- 


B{p-1) - A{0 -1) 

1 + B{p-1)-A{0-1) 


)■ 


(2.34) 


= ^l_A{0 + 0r)+Bp + p^’^{p,0))k^ 


A,B 


where 


{0 + Or) d^,cl>^’^{p, 0) = cf'^(p, 0)) fcf, 

pf’^(p, 0) = -Ab({p - l)i±^ + + A^{0 - 1) 

V XA,B XA.B-' 

B 


kf’^ = 


XA,B 

( 

B 

1 


0 - 1 ' 

XA,B 

B 


1 -{- Or 
XA,B 

B 


B 


2 P 


- 1 


XA.B 


, , u 

3. , „ 13/12/ “ ^4 (1 13/12 

^A.B ^A.B 


^ Xa.B 

ABi -XA,BP{0 + 0r){2 + B{p-2)-2A{0-l)) 

c, ’ (p, 0 ) = -^ 

{1 + B{p-1)-A{0- 1)) (1 -B-A{0- 1)) 

k^'^ = (kt'^)-\ 
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Note that due to (j2.32|l we have B 0 it A, B tend to zero. Thus, our 

example fulfills the requirements ()2.26|l - (j2.28ll as well as the requirements for (I2.22|l 
and (I2.23|l . 


2.2. Constitutive limit. Now we want to show that we are able to obtain the 
Oberbeck-Boussinesq approximation as a constitutive limit of (12.251) as A, B tend 
to zero. We assume that the requirements (I2.26I) ~ (I2.28I) as well as the requirements 
for (12.221) and (12.231) are fulfilled. We consider solutions v = 9 = 9^'^ and 

p = of (j2.25p in the following sense: 

Let n C be a given sufficiently smooth bounded domain and / = (0, T), T > 0, 
be a given time intervall. We set Q = I x il and assume that V/ € L^{Q) and 
appropriate boundary and initial conditions for v and 9 are given. The system 
(I2.25|) is satisfied in the sense of distributions over Q and v, 9 and p are uniformly 
bounded in D L°°{Q) with respect to A, B. Finally, we assume that 

B = o{A). 

From our assumptions follows that there exists v, 9,p £ such that 


- V 

weakly in 

9 A.B Q 

weakly in W^'‘^{Q), 

pA,B^p 

weakly in W^’^{Q) 

^ V 

strongly in L‘^{Q), 

9A,b Q 

strongly in L'^{Q), 

p^’^ -A p 

strongly in L'^{Q). 


These convergences and the assumptions ()2.26p - (j2.28l) as well as the requirements 
for (12.221) and (|2.23l) immediately imply that the limit as A, B tend to zero in 
(12.251) 1 -^ yield for all £ C^{Q) 


'Q 


'Q 


Co dt9 if) A cqv ■ V9 Ip £- 


PrRe, 


divv ■ (f dx dt = 0, 
■ • Vi/) dx dt = 0. 


(2.35) 


To treat (j2.25|) o we add on both sides and obtain for all ip £ C^{Q) 


I 


- A (9^’^ + 9r) + Bp^'^ + 


X ■ ipdxdt 

——■ Da? H-div v^’'® div p> dx dt 

Reu ReA 


9^’^))x 


- j -^ ’^/(l - ^r) - -^f) div if dx dt 

J Q 
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From this we deduce for A, B tending to zero for all if S C^{Q) with divy? = 0 


[ (9tv + [Vv]v) ■(pdxdt+ [ 
JQ Jq 

= [ {l — f ■ ipdxdt. 

Jo 


Dv ■ Dip dx dt 


(2.36) 


The system (I2.35I1 . (12.3611 is the weak formulation of the celebrated Oberbeck- 
Boussinesq approximation. 


2.3. Formal expansion. Our approach enables us also to recover the results in [12] 
and |5] for dissipation number Di = 0. To this end we assume that the parameters 
A, B are fixed but sufficiently small. Moreover, we assume that A is smaller than B 
in the sense that A^ < B < A. This requirement is fulfilled in typical applications. 
In this situation we formally expand the non-dimensional quantities v, 9 and p in 
(I2.25|) into power series with respect to the perturbation parameters A, B in the 
form 

OO OO CXD 

V = ^ A^B’^Wm^k, 6 = ^ A^BHm^k. P=Y. 

m^k—0 m^k—0 7n,k—0 

We attach the boundary conditions of the quantities to the zero order term and set 
the boundary conditions to be zero for the higher order terms. Similarly we expand 
the quantities co,kf’^, z = 1,2, and functions af’^{p,9), pf’^{p,9), cf’^{p,9), 
p^’^{p,9), z = 1,2, where we also insert the series for p and 9 in the arguments. 
We replace all these quantities in the system (I2.25|) by the corresponding power 
series. At level (0,0) we see from (12.2511 thaid 


co^o.o “ 


1 


Pr Rou 


divvo,o = 0, 
7Vpo,o = V/, 

A0O.O = 0. 


At level (1,0) we see from (12.251) 9 that 

vo,o - T^divDvo.o + qVpi.o = -(6*0,0 - 6*r)V/ -f kippVf. 

Keu 


Thus, setting v := vq^q, ^ := 6*o,o andp := 7Po,o + A(7Pi,o —6*r/ —fci;i,o/) we obtain 


divv = 0, 


a(v — ——divDv) -f Vp = (1 — A9)Vf, 
\ Re,, / 


co9 - 


1 


Pr Re,, 


-A9 = 0. 


(2.37) 


This system is the same as the one obtained in m and 13 for Di = 0. The approx¬ 
imation combines different levels in the temperature and the velocity equation and 
the velocity has zero divergence. 


^To simplify the notation we use the superposed dot only at the level (0, 0). 
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